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We use the linear a model with two flavors of quarks to study the phase diagram at finite tem- 
perature and baryon chemical potential as a function of the vacuum pion mass. Our calculations 
include thermal fluctuations of both the bosonic and fermionic fields. An interesting phase struc- 
ture occurs that results in two, one, or no critical points depending on the value of the pion mass. 
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The physical pion mass is small but not zero. In consequence, the conventional wisdom is that 
for two flavors of quarks there is no true thermodynamic chiral phase transition at finite temperature 
T and zero baryon chemical potential \i. Instead, there is expected to be a curve of first-order phase 
transition in the \i-T plane that terminates in a second-order phase transition at some critical point 
(fi c ,T c ). The location of the critical point obviously depends on the physical (vacuum) pion mass. 
This topic has been under intense theoretical study using various effective field theory models, 
such as the Namu Jona-Lasinio model [|l|, ||, ^J, a composite operator model [ffl|, a random matrix 
model |Bj|, a linear a model [Eh, an effective potential model and a hadronic bootstrap model 
[^], as well as various implementations of lattice QCD [|| ^, [T(| [Tlj ]. Reviews of the subject were 



presented in the past few years [12]. It is also of great interest because collisions between heavy 
nuclei at medium to high energy, such as at the future Facility for Antiproton and Ion Research 
(FAIR), or possible low energy runs at the Relativistic Heavy Ion Collider (RHIC), may provide 
experimental information on the phase diagram in the vicinity of a critical point. 

In this work we study the phase diagram of the linear a model coupled to two flavors of 
identical mass quarks. Whereas the a model is an oft-used effective model that represents some of 
the essential features of the chiral dynamics of QCD, the reason to couple the fields to quarks is less 
obvious. One argument is based on the existence of the critical point itself. If a critical point exists, 
then one can go around it without crossing the curve of the first-order phase transition. The effective 
degrees of freedom should not change too much in following such a path. Therefore, if constituent 
quarks are considered to be reasonably useful degrees of freedom on the higher temperature side 
then they should be useful on the lower temperature side too. For this reason we use quarks to 
carry the baryon number but acknowledge the resulting uncertainty in the results. The present 



work extends that of Refs. [g, |13L |14[ |15|, |16|] in several ways. We include thermal fluctuations of 
the meson and fermion fields, which can be important at finite temperature when the magnitude 
of the fluctuations becomes comparable to or greater than the mean values. We scan vacuum pion 
masses from zero to over 300 MeV, which may be particularly useful for comparison with lattice 
gauge theory with different quark masses. Even when all other parameters of the model are fixed, 
we find an interesting phase diagram that may have two, one or no critical points depending on the 
value of the vacuum pion mass. 
The Lagrangian is 

& = I (d^) 2 + I (d^f ~ U(cr, 7r)+i>[i^-g(a + i 75 r • tt)] V , (1) 



where 



U(a,ir) = -(a 2 + ir 2 -f 2 f-Ha (2) 



in an obvious notation. The SU(2) l xSU(2)jj chiral symmetry is explicitly broken by the term Ha, 
which gives the pion a mass. The parameters in the Lagrangian are constrained by fixing the pion 
decay constant /„. = 92.4 MeV, the a mass m a = 700 MeV, and the quark mass is set to one-third 
of the nucleon mass m q = 313 MeV. The vacuum pion mass is varied from to m a /2 = 350 
MeV. To create an effective mesonic model, we integrate out the quark degrees of freedom in the 
usual way [|T7j], such that 

In iT quark = In detD, (3) 
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where D is the inverse quark propagator. This leads to an effective Lagrangian 
where 



^eff = ^ {d^Y + - (d^a) 2 - U eS (a, tt) , (4) 



T 

Ues(a, tt) = U (a, tt) - — In 5° quark (cr, tt) (5) 

is the effective potential. We decompose the scalar field into a condensate v plus a fluctuation A. 
The thermodynamic potential is computed from 

O = (C/eff) - \m 2 a (A 2 ) - ^m 2 (tt 2 ) + ^ + Or , (6) 

where £2^ and Q.^ are the independent particle contributions from the a and pion quasiparticles. 
Angular brackets refer to thermal averaging. The condensate is determined self-consistently by 



while the quasiparticle masses are determined self-consistently by 

2 _ / d 2 U e ff \ 2 _ / ® 



(7) 



dA 2 / 77 \ d 



< = ( "^P ) < = { -^r-r ) • (8) 



77T 



The thermal fluctuations are determined by 



o 



and 

.2 



„ aar 3 f p 2 1 
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The techniques for evaluating the thermal average (J7 e ff) and its derivatives are nontrivial and were 



developed in Refs. [13, 14, 15]. The resulting equations must be solved (very accurately) using 
numerical methods. All thermodynamic identities were verified both analytically and numerically. 

An example of the pressure versus chemical potential at fixed temperature is shown in the left 
panel of Fig. 1 . For the higher temperature of 80 MeV there is only one self-consistent solution. 
For the lower temperature of 50 MeV, there is a unique solution at large \i and another unique 
solution at small ji. For a range of \i centered about 900 MeV there are three solutions: One is 
associated with a continuation of the low-density phase, one is associated with a continuation with 
the high-density phase, and the third solution (not shown in the figure) is an unstable phase. The 
point where the two curves cross is the location of the phase transition. In this example it is first- 
order since there is a discontinuity in the slope dP{ji,T) / 'dfi = ub- Where each curve terminates 
is the limit of metastability for that phase. The thermodynamically favored phase is the one with the 
largest pressure. At some temperature between 50 and 80 MeV the slopes are equal at the crossing 
point, there are no metastable phases, and the second derivative is discontinuous. This is indicative 
of a second-order phase transition. The location of this point in the /j,-T plane is the critical point. 
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The phase diagram for a sampling of vacuum pion masses is shown in the right panel of Fig. 
1. Generically there is a curve of first-order phase transition starting on the fi axis and arching to 
the left. This curve terminates at a critical point. For m^^c = 321 MeV the critical point sits on 
the T = axis, and for m IiVac > 321 MeV there is no phase transition at all. Of course the precise 
numbers depend on the constants in this model, such as the vacuum a mass and what value one 
assigns to the coupling of the quark field to the a field, but the results are in line with expectations 

& 

Something very interesting happens as the vacuum pion mass is decreased from 50 MeV to 
0. The left panel of Fig. 2 shows the phase diagram for a vacuum pion mass of 35 MeV. There 
are now two critical points! There is a line of first-order phase transition beginning on the /j, axis 
and arching to the left to end at a critical point of \x c \ 725 MeV and T c \ m 92 MeV and there is 
another line of first-order phase transition beginning on the T axis and arching to the right to end at 
a critical point of fj, c 2 ~ 240 MeV and T C 2 « 137 MeV. The latent heat curves are shown in the right 
panel of Fig. 2. In the vicinity of \i = 575 MeV and T = 1 10 MeV the latent heat gets pinched to 
zero for some critical value of the vacuum pion mass between and 35 MeV. 

The results shown in the right panel of Fig. 1 illustrate the conventional view on the nature 
and location of the critical point. The result shown in the left panel of Fig. 2 is unconventional or 
exotic. To the best of our knowledge, this is the first time a model calculation has resulted in such a 
phase diagram. It is well known that the nature of the phase transition or crossover for two flavors 
of quarks is sensitive to such details as the strength of the axial U(l) anomaly and the value of the 



vacuum a mass [ ]1 8| , |19| , |20| ]. Even for two flavors of massless quarks with \i = the order of the 
transition is not known. In particular, the authors of Ref. [21] found that their Monte Carlo results 
were substantially consistent with a first-order transition with respect to scaling of the specific heat 
and the chiral condensate, but not with respect to the chiral susceptibility. It cannot be expected that 
a model as simple as the one studied here can make definitive predictions for what actually happens 
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Figure 2: Left panel: Phase diagram for a vacuum pion mass of 35 MeV. There are two critical points. 
Right panel: Latent heat along the coexistence curve for illustrative values of the vacuum pion mass. Note 
the second region in the lower right corner for the 35 MeV mass. 



in QCD. However, it can serve to illustrate the possibilities. One such possibility is sketched in Fig. 
3. There are critical surfaces which delineate regions of first order phase transitions. Some lattice 
calculations JToj ] show the surface bending towards the /i-axis. This cannot continue indefinitely, 
and it is conceivable that the surface will eventually curve away from it. This would indicate two 
critical points for certain ranges of the quark masses. 



Real world 
Heavy quarks 
N,=3 




Figure 3: Sketch of what might happen with three flavors of quarks. Due to a change in the curvature of the 
surface more than one critical point might arise. Figure courtesy of Ph. de Forcrand. 
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Further work with this and related models can be easily identified. Changing the other param- 
eters in the model, such as the vacuum a mass or the strength of the coupling of the quarks to the 
chiral fields, will obviously change the particular values of the pion masses for which two, one or 
no critical points exist. Recent studies in the Nambu Jona-Lasinio model with varying vector cou- 
pling strength show analogous critical behavior [22]. It would be interesting to extend the a model 
to include strangeness and to vary both the vacuum pion and kaon mass using the self-consistent 
techniques of the present paper. It should be straightforward to calculate shear and bulk viscosities 
and thermal conductivity along the lines of Refs. [23, 24]. These transport coefficients could be 
used in fluid dynamic modeling of heavy-ion collisions at FAIR and for low-energy runs at RHIC. 
More details on this work have been, or will be, presented elsewhere 26]. 
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